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Introduction 

It  is  the  puTi'Ose  of  this  report  to  give  a  detailed  account  of  the  elements 
of  the  theory  of  Hill's  equation,   to  review  briefly  the  more  sophisticated  recent 
results  and  to  provide  a  guide  to  the  extensive  literature  on  special  cases  and 
their  applications  to  mechanical  and  electrical  problems. 

The  first  part  of  this  report  deals  with  the  general  theory  of  Hill's  equation. 
In  Chapter  1  and  in  the  first  three  sections  of  Chapter  2  the  elementary  parts  of 
the  theory  are  developed  and  complete  proofs  of  all  results  are  given.     Prerequisites 
for  these  sections  are  the  elementary  theory  of  ordinary  linear  differential  equations 
and  standard  results  in  the  theory  of  functions  of  a  complex     variable.     The  last 
three  sections  of  Chapter  2  have  the  character  of  a  review  and  proofs  of  the  quoted 
results  have  been  omitted. 

A  brief  historical  introduction  to  Hill's  equation  can  be  found  in  Moulton,   1930. 

The  references  in  the  present  report  do  not  always  give  the  earliest  source 
for  a  theorem.     In  particular,    certain  results  usually  referred  to  as   ''Haupt's 
theorem* •   have  been  quoted  under  this  name  although  they  were  first  derived  by 
Liapounoff,   190?. 

For  some  interesting  approaches  to  the  theory  of  Hill's  equation  which,    so  far, 
have  been  applied  to  a  special  case  (Mathieu's  equation),    see  Chapter  20  in  Whittaker- 
Watson,   1927. 
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1,  Basic  Concepts 

1«1«     Preliminary  remarks 

Any  homogeiieous  linear  differential  equation  of  second  order  with  real 
periodic  coefficients  can  be  reduced  to  an  equation  of  Hill's  type«     The  specific 
question  which  arises  in  the  theoi^r  of  Hill's  equation  is  the  problem  of  the 
existence  of  periodic  solutions*     This  problem  has  many  features  in  common  with  the 
ordinary  Sturm-Liouvilie  problems,   and  in  certain  cases  it  can^in  fact^be  reduced 
to  ordinary  boundary  value  problems  of  the  Stura-Liouville  type,    (see  Section  1.3 )♦ 
However,   in  general  such  a  reduction  is  not  possible,   and  imposing  the  periodicity 
requirements  on  a  solution  of  the  differential  equation  leads     to  phenomena  different 
from  those  resulting  from  the  imposition  of  a  homogeneous  boundary  condition  of  the 
Sturm-Liouville  type.     Thus,   the  differential  equation  can  have  two  linearly 
independent  periodic  solutions  but  it  cannot  have  two  linearly  independent  solutions 
satisfying  the  same  homogeneous  boundary  conditions.     Futhermore,  the  value  of  the 
period  of  the  solution  (which  is  a  multiple  of  the  period  p  of  the  coefficients) 
plays  an  important  role  in  the  discussion  of  periodic  solutions.     In  a  certain 
sense,   only  the  solutions  of  period  p  and  2  p  are  of  interest.     We  shall  now  proceed 
^with  a  detailed  presentation  of  some  basic  theorems  and  their  proofs. 

As  references  to  the  general  theory  of  Sturm-Liouvilie  (self -adjoint)  boundary- 
value  problems  we  mention  Courant  and  Hilbert,  1953  and  Coddington  and  Levinson,   1955« 


1.2     Floqiiet's  Theorem 

Let  Q(x)  be  a  (real  or  complex  valued)  function  of  the  real  vairlable  x  defined 
for  all  values  of  x.     We  assume  that  Q(x)  is  piecewise  continuous  in  every  finite 
interval  and  that  it  is  periodic  with  minimal  period  n.     By  this  we  mean  that  for 
ail  X 

(1.1)  Q(x+n)  -  Q(x), 

and,  that  if  p  is  a  number  with  0  <  p  <  n,  then  there  exists  at  least  one  real  intervai 
1  such  that  Q(x+p)  /  Q(x)  for  x  £  I. 

If  Q(x)  has  the  properties  stated  above,  then  the  differential  equation 

(1.2)  y"  +  Q(x)y  -  0 

has  two  continuously  diff  erentiabie  solutions  y,  (x)  and  yp(x)  which  are  uniqjiely 
determined  by  the  conditions: 


-  3  - 

y^(0)  -  1,  y^(0)  -  0,   y^CO)  -  0,   y^(0)  -  1. 

These  solutions  are  referred  to  as  normalized  solutions  of  (1«2). 

Before  stating  Floquet 's  theorem  we  must  define  the  notions  of  characteristic 
equation  and  characteristic  exponent  associated  with  (1.2 ).  Thus,  the  characteristic 
equation  is  the  equation 

(1.3)  9'  -  |yi(n)  +  y^(it)]  p  +  1  -  0 

and  the  characteristic  exponent  a  is  a  number  which  satisfies  the  equations 

(l.U)  9xp  ian  ■  P]^  »  eocp(-ian)  ■  P2  » 

where  p^  and  p.  are  the  roots  of  the  characteristic  equation  (1*3 )« 

It  is  clear  that  a  is  defined  up  to  an  integral  multiple  of  2.  Also 

2  cos  an  -  y^(n)  +  y^(n).  Finally  p^p^  ■  1. 
We  can  now  state  Floquet 's  Theorem: 

1  .  If  the  roots  p.  and  p^  of  the  characteristic  equation  (1»3)  M'e  different  ' 
from  each  other,  then  Hill's  equation(l»2)  has  two  linearly  independent  solutions 

f^(x)  -  e^^*p^(x),  f^Cx)  -  e'^^'^p^Cx), 

^ere  p,(x)  and  Po(x)  are  periodic  with  period  n. 

2  .     i£.  Pi   ■  Pp»   then  equation  (1.2)  has  a  non-trivial  solution  which  is 
periodic  with  period  n   (vh  en  p.   -  p_  ■  1)  or  2n   (when  p,   -  p-  f  -1),     Let 
p(x)  denote  such  a  periodic  solution  and  let  y(x)  be  another  solution  linearly 
independent  of  p(x)»     Then 

y(x+n)  -  p     y(x)   ♦  9  p(x),         9  constant, 

and  Q  »  0  is  equivalent  to 

y^C")  +  /^(n)  "12,     y^Cn)  -  0,     yj^(n)  -  0. 

Before  starting  with  the  proof  of  Floquet 's  theorem  it  may  be  appropriate  to 
discuss  its  significance.  Thus,  let  p^  f  pp.  If  a  is  real,  then  there  exists  an 
upper  bound  Fi   for  the  absolute  value   |y(x)|   of  every  solution  of   (1.2)  and  M  depends 
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on  the  initial  conditions  for  y  and  not  on  x.     If  a  is  not  real,   then  there  exists  a 
non-triTlal  unbounded  solution  y(x)  of   (1»2).     If  p-   ■  p^,  then  for  all  solutions 
of   (1.2)  to  he  bounded  it  is  necessary  and  sufficient  that 

y^M   +  y^(")  "12,     j^in)  -  0,     jj^in)  -  0. 

Whenever  all  solutions  of   (1.2)  are  bounded  we  say  that  they  are  stable.     Other- 
wise we  say  that  they  are  unstable. 

The  solutions  of  period  n  and  2n  play  an  exceptional  role  as  is  seen  from  the 
following 

Corollary  to  Floguet's  Theorem.     If (1.2)  has  a  periodic  non-trivial  solution 
with  period  nn,  n  >  2,.  but  no  solution  with  period  n  or  2n,   thwi  all  solutions 
are  periodic  with  period  roi. 
Indeed,   our  assxunption  implies  that  p.   ^  p.  so  that  every  solution  y  of   (1.2)  is 
of  the  form 

y  -  iif^M  *  vf^ix). 

If  one  such  solution  is  periodic  with  period  nn,   then  y(x+nH)  ■  iicf^    ♦  vcf^  -  y(x) 

where  c  -  exp(iann),   c  »  exp(-iann).     Since  f     and  f^  are  linearly  Independent, 

c  ■  c  "  1,     Therefore,   na  is  an  even  integer,   and  both  f^  and  f.  are  periodic  with 

period  nn. 

Proof  of  Floquet's  Theorem.  If  y(x)  is  a  solution  of  (1*2 ),  then,  obviously,  y(x-hx) 

is  also  a  solution  of  (1.2).  In  particular,  y,(x+n)  and  y-(x+n)  are  solutions  of  (1.2). 

Since  y,(x)  and  y^M   form  a  basis  for  the  set  of  all  solutions  of  (1.2),  it  must  be 
possible  to  express  y, (x+n)  and  y  (x+n)  as  linear  combinations  of  y,(x)  and  yp(x).  We 
find  easily  that 

y^(x+Ti)  -  y^(n)y^(x)  +  y^(n)y2(x) 


(1.5) 


( 


y2(x+n)  -  y2(n)y^(x)  ♦  y^(n)y2(x) 


Assume  now  that  y(x)  f  0  is  a  solution  of  (1.2)  such  that 

(1.6)  y(x+n)  -  py(x) 

for  some  constant  p.  If  y(x)  ■  c-y^(x)  +  C2yp(x),  then  it  follows  from  (1«6)  that 
c,  and  Cp  must  satisfy  the  system  of  linear  equations 


(1.7) 
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y|(n)c^  +   (y^(n)  '1^)0^  -  0 


Conversely,   if   (1»7)  is  satisfied,   y(x)  satisfies   (1.6),     Now,  the  necessary  and 
sijfficient  condition  for   (!•?)  to  have  a  solution  c,,   c-  such  that  not  both  c-   and 
Cp  vanish  is 


(1.8) 


yi(")-p  >  y2(") 


y{W  ,  y^(n)-p 


0  . 


Since,   for  all  x,   the  Wronskian 


yi(x)y^(x)  -  y2(x)y«i(x)   «  l, 

equation  (1»8)  is  identical  with  tl-ie  characteristic  equation   (1»?).     Thus,   if  p  ■  p^ 
is  a  root  of   (1.8),  we  can  find  c^   and  c^  rnich  that  y  «  c,y^   ♦  ^^y?  ^  ^  ^^   such  that 
y  satisfies   (1.6).     Obviously,   if  (1.6)  is  satisfied,  we  may  write 

y  -  y(x)  -  exp(iax)  p^(x)  -  f^(x) 

where  exp(ian)  ■  p^  and  where  p,(x)  is  a  periodic  function  of  x  with  period  n.  Suppose 
now  that  (1»8 )  has  a  second  solution  p  ■  p^  /  p. .  We  may  use  p.  for  the  construction 
of  a  solution  y  -  f2(x)  ^0  of  (1.3)  such  that  f2(x+n)  -  p2f2(x).  We  observe  that  f 
and  fp  are  linearly  independent.  Otherwise,  we  could  find  constants  X^  and  X  not  both 
of  which  vanish  such  that 


But  then  we  would  have 


\f^(x)  ♦  X2f2(x)  B  0. 


X^f^(x+n)  +  X2f2(x+n)  -  >^Pii'i(x)   +  X2P2i'2(x)  =  0» 


Since  \f-,   and  ^fp  do  no*  both  vanish  identically,   the  last  two  equations  ai  e 
compatible  only  if  p-   ■  p.,   which  we  have  excluded.     This  proves  Floquet's  theoran  in 
the  case  where  p,   /  p^,. 

Since  p^pg  -  1,    either   |p^|   ■   Ipgl   ■  1  or  at  least  one  of  the  numbers   |p.|,    |p-| 
exceeds  1,     In  the  first  case,   we  have  stability^in  the  second  case  instability  of  the 
solutions  of   (1.2),   provided  that  p^   T^  Pr,* 
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If  P-.  ■  p«,  we  still  can  construct  one  solution  y^  (x)  of  (1«2)  such  that 

Since  p^  ■  p«  and  p^p-  ■  i  implies  that  p^  "  1  1>  y-i  is  obviously  periodic  with 
period  n  or  2n.  In  order  to  find  the  properties  of  a  solution  y^Cx)  which  is  linearly 
independent  of  y-,  assume  first  that  y2(n)  /  0.  Then  we  may  choose  (cf.  (1«7)  and  (l.JJ 


y*(x)  -  7^{n)y^ix)   ♦  lpi-7i(«)j  72 U) 


72 (x)  -  y2(x) 
and  we  find  from  2p^  -  y^(n)  +  y^(n)  that 

y*(x+n)  -  pj^ygCx)  +  yj^(x). 
Similarly,  if  yp(n)  "0,  we  may  choose 

y^Cx)  •  ygCx),  ygCx)  -  y^^Cx). 

Since  y^(n)y^(n)  -  y^{n)y^in)   ■  1  it  follows  from  y^in)   -  0  and  y^in)   ♦  y^(n)  •  Zp^ 
that  y, (n)  ■  yA(n)  "  p,  and  therefore  we  have  from  (1.5)  that 

y*(x+n)  -  p^y^(x) 

y*(x+it)  -  p^ygCx)  ♦  y^(ii)yj^(x). 

This  proves  Floquet's  theorem  in  all  details. 

As  a    rather  obvious  consequence  of  Floquet's  theorem  we  mention  the  following 
Stability  Test.     The  solutions  of   (1.2)  are  stable  if  and  only  if 
7i(n)   ♦  7j(i»)  is  real  and 

ly^Cn)  ♦  7^(1)1  <  2 


or 


and 


y^Cn)  ♦  y|(")  -  1  2 


ygCn)  ■  y{(")  ■  0' 


Proof.     If  p^  f  P-,   then  stability  is  equivalent  to  a  ;*  0,   a  real,   which^in  turn,  is 
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equivalent  to  y, (n)  +  y^(n)  being  real  and  in  absolute  value  <  2.  If  p-  ■  pp,  then 
stability  is  equivalent  to  y, (n)  +  y^(n)  "  ♦  2  and  yp(n)  -  yJ (n)  ■  0. 

1.3  The  symmetric  case  Q(x)  ■  Q("x) 

If  in  (1»2)  the  function  Q(x)  is  even,  i.e.,  if 

(1.9)  Q(x)  -  Q(-x) 

it  is  possible  to  establish  relations  between  the  values  of  y^ ,  y.,  yJ,  y »  at  x  »  n/2 
and  at  X  ■  n  and  these  relations  allow  a  more  detailed  classification  of  the  solutions 
of  period  n  and  2n.  We  summarize  our  results  by  stating 

Theorem  1,1,     Let  y, (x)  and  y^Cx)  be  the  normalized  solutions  of  (1»2)  and 
assume  that  Q(x)  satisfies  (1«9)«  Then  the  following  relations  holdt 

(1.10)  y^(n)  -  2y^(J)y'(J)  -  1  -  1  ♦  2y|(J)y2(5-) 

(1.11)  y2(Ti)  -  2y^{^)Y^q) 

(1.12)  y'(n)  -  2y^(J)yj(|) 

(1.13)  y^(n)  -  y^(n). 

Theoran  1.2.     If  the  conditions   of  Theorem  1.1  are  satisfied,   then  there 
exists  a  non-trtvial  periodic  solution  of   (1.2)  vriiich  is 

(1)  even  and  of  period  n  if  and  only  if  y-J  (^)  ■  0 

(2)  odd  and  of  period  n  if  and  only  if  jji)  •  0 

(3 )  even  and  of  period  2n  if  and  only  if  y. (i)  ■  0 

(U)  odd  and  of  period  2n  if  and  only  if  yMw-)  "  0. 

Periodic  solutions  of  period  w  or  2n  are  necessarily  multiples  of  the 
normalized  solutions  y, (x)   and  y^{x). 

Proof  of  Theoron  1.1.     If  Q(x)  is  even  and  if  y(x)  is  a  solution  of  (1.2),   then  y(-x) 
is  also  a  solution.     Since  the  initial  conditions  for  y,(-x)  and  y, (x)  coincide  and, 
similarly,  those  for  yAx)  and  -yp(-x)  are  identical,   it  follows  that  y,  (x)  is  even 
and  y2(x)  is  odd.     Therefore  we  find  from   (1.5)  for  x  -  -ti/2; 

(l.lU)  7^{^)  -  yi(n)y^(J)  -  j^My^q) 
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(1*15)  ygCj)  -  y2(")yi(5')  -  y^(n)y2(j). 

Obviously,  yJ(x)  is  odd  and  yi(x)  is  even.  Using  this  fact,  we  find  from  (1.5) 
by  differentiating  both  sides  with  respect  to  x  and  by  putting  x  ■  -n/2  that 

(1.16)  y^q)   -  -yi(n)y'(J)  +  y'(n)y^(J) 

(1.17)  y|(J)  •  -y2(n)yi(j)  ♦  y^('t)y^(J)  . 

We  may  treat  (l»lU)  to  (lol7)  as  a  system  of  linear  equations  for  y, (n),  yp(n), 
y-J(")»  y^W*  Solving  the  equations  for  these  quantities,  we  arrive  at  Tnecrem 
1.1  if  we  observe  that 

y2^(J)y^(5')  -  y^(^)y2(f)  "  i* 

Proof  of  Theorem  1»2.     It  is  easily  seen  that  an  even  solution  of  (1.2 )  is  a  multiple 
of  y,  (x)  and  thiat  an  odd  solution  must  be  a  multiple  of  yp(x)» 

To  prove(l)  let  us  assxime  that  y(x)  is  a  non-trivial,  even,  periodic  solution 
of  (1*2)  with  period  n.     Then  y, (x)  is  also  peidodic  wLth  period  n  and  the  same  is 
true  about  yJ(x).     Thus  yiij)  ■  y{(-  x).     Since  y4(x)  is  an  odd  function  it  follows 
that  yj'(|)  -  -yji(-  J).     Consequently,  yjf(j.)  -  0.     Conversely,  if  yj^(J)  -  0,  then 
y»(-  *•)  ■  0.     Also,   y-,  (-  ^)  ■  yi(^).     Therefore  y,  (x)  satisfies  the  same  conditions 
at  X  »  -  J  and  at  x  -  »  .     From  this  and  from  the  periodicity  of  Q(x)  it  follows 
that  y,  (x)  is  periodic  with  peii-od  n. 

The  proof  of  (2)  is  entirely  analogous  to  the  proof  of  (l). 

In  proving   (3)  we  show  that  the  values  of  y,  (x)  and  yJ(x)  at  x  ■  -  i-    differ  in 
sign  from  the  values  of  y-i  (x)  and  y^(x)  at  x  «  ^  .        This  shows  that  y^(x*r.)  - 
-3''-i(x)  and  therefore  y,(x+2n)  ■  y,(x). 

The  proof  of   (U)  is  analogous  to  the  proof  of   (3)» 

2«  Characteristic  Values  and  Discriminant 

2.1.  Characteristic  values  and  intervals  of  stability 

In  this  and  in  the  following  sections  we  shall  study  Hill's  equation  in  its 
standard  form 

(2.1)  y"  +  jx  ♦  Q(x)T  y  -  0, 

where  X  is  a  parameter  and  where  Q(x)  is  a  real  peiriodic  function  of  x  with  period  n. 
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Unless  otherwise  stated  we  shall  assume  that  Q(x)  is  two  times  differentiabie 
for  all  X.  T'liis  is  a  rather  strong  assumption,  and  occasionally  we  shall  treat 
examples  where  QCx)  is  not  even  continuous. 

Let  y^  and  y-  be  the  two  linearly  independent  solutions  of  (2.1)  which  we 
defined  by  simple  initial  conditions  in  Section  1.2.  To  emphasize  their  dependence 
on  X  we  shall  sometimes  write  y^(x,X)  and  ypCxjX)  instead  of  y^(x)  and  ygCx).  One 
of  our  main  problams  w:.ll  be  the  determination  of  those  values  of  X  for  which  the 
solutions  of  (2.1)  are  stable.  This  problem  is  idaitical  with  the  problen  of 
determining  those  values  of  X  for  which  equation  (2.1)  has  a  solution  of  period  n 
or  2n. 

The  following  theoran  due  to  0.  Haupt  (191ii,  1918)  connects  the  two  problems. 
Theorem  2.1.  (Oscillation  Theorem)  To  every  differential  equation  (2.1), 
there  belong  two  monotonically  increasing  infinite  sequences  of  real  ntimbers 

(2.2)  X^,  X^,  X^,  ••. 
and 

(2.3)  X|,  X^,  Xj,  X^,  ... 

Such  that  (2.1)  has  a  solution  of  period  n  if  and  only  if  X  ■  X  ,  n  ■  0, 
1,  2,  ...  and  a  solution  of  peirLod  2n  if  and  only  if  X  •  X»,  n  -  1,  2. 
.3 , . . .  .  The  X  ,  X »  satisfy  the  inequalities 

(2.3)         X^<X|<X^<X^<X2<Xj<Xj^<X^<X^<,,, 

and  the  relations 

lim       X"-"-  -  0     ,  lim         (X')"''-  -  0. 

Ac)  n— >oo  n-E>oo 

The  solutions  of   (2.1)  are  stable  in  the  intervals 

(2.6)  (X^,    X^),    (X^,    X^),    (Xg,    Xj),    (Xj^,   X^),    ...      . 

At  the  endpoints  of  these  intervals  the  solutions  of  (2.1)  are,  in  general, 
unstable.  This  is  always  true  for  X  ■  X  •  The  solutions  of  (2.1)  are  stable 
for  X  -  X   -  2£  ^  ■  ^  +2  if  and  only  if  X-  ^  «  X-  p*  and  they  are  stable 
£21  ^  ■  Mn*l  SL^'^iz   if  and  only  if  X-^^^  -  X'^^^. 

For  complex  values  of  X  (2.1)  has  always  unstable  solutions. 
In  order  to  be  able  to  refer  briefly  to  the  assertions  of  Theorem  2.1,  we  shall 
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use  the  following  definitions: 

The  real  numbers  X  shall  be  called  characteristic  values  of  the  first  kind 

n  '  ' — ——————— 

of   (2.1)  and  the  X»  shall  be  called  characteristic  values  of  the  second  kind.     The 

n  ,.«,,^__.__.«-^____--.— -—^-— —————— - 

intervals   (2.6)  on  the  real  X-axJ s  shall  be  called  intervals  of  stabilityj  an  end 

point  of  such  an  interval  shall  belong  to  it  if  and  only  if  (2.1)  has  stable  solutions 

for  the  corresponding  value  of  X,     Similarly,  we  shall  talk  about  intervals  of 

instability.     Both  the  intervals  of  stability  and  of  instability  are  ordered  in  a 

natural  manner.     We  shall  always  disregard  the  interval  of  instability   (-oo,   X^). 

Therefore,  the  first  interval  of  instability  will,  in  general,  be  the  interval 

(X»,    X*).     Observe  that,  according  to  Theorem  2.1,   neither  an  interval  of  stability 

nor  an  interval  of  instability  can  ever  shrink  to  a  point.     The  intervals  of  stability 

can  never  disappear,  but  two  of  them  can  combine  to  a  single  one  if  X-     ,  ■  ^-+2  '^^ 

\X     ,  -  X»     -.     However,  the  intervals  of  instability  may  disappear  altogether  (e.g., 
^n+1         2n+2 

if  Q(x)  is  a  constant).     It  is  a  deep  result  due  to  Q.  Borg,  19U6,   that  for  a  non-con- 
stant Q(x)  there  exists  at  least  one  interval  of  instability.     As  we  shall  see 
later  (Section  2.U)  there  may  exist  only  one  interval  of  instability. 
Proof  of  Theorem  2.1.     We  shall  exclude  first  the  possibility  of  stable  solutions 
of  (2.1)  in  the  case  of  a  complex  value  of  X.     Assume  that  X  ■  tx  +  iv,  where  n,  v 
are  real  and  v  /  0.     Let  y  -  u  +  iv  be  a  solution  of  (2.1)  which  is  of  the  type 

(2.7)  y(x)  «  e^"""  p(x)  -  u  +  iv 

where  o  is  real  and  where  p(x)  is  periodic  with  period  n.     According  to  Floquet's 
Theorem,  such  a  solution  y(x)  exists  if  we  have  stability  for  the  solutions  of  (2«1). 
3y  splitting  (2.1)  into  its  real  and  imaginary  parts,  we  find 

u"  +  L  ♦  Q(x)     u  •  vv 

(2.8) 

V  •  +     p.  +  Q(x)     V  ■  -vu     . 

If  w«  multiply  the  first  of  the  equations  (2.3)  by  v  and  the  second  by  u  and  fonn 
the  difference  of  the  resulting  equations  we  find  that 

2    2 
ut»  V  -  v'  'u  ■  V  (u  +  V  )  , 

or,*  upon  integrating,   that 


(2.9)  u'v  -  v'u  -  V  f     |u^(t)  +  v^(t)j  dt  ♦ 

/  o 


^(t)  ♦ 

v^(t)     dt   . 

e  u     ♦ 

V    ■    |p       and 

th 

erefore,   for  n  « 

1, 

2, 

3, 

0  •  •  ) 

u2(t) 

♦  v^(t)     dt  - 

n 

n 
1      lp(t)  2  dt. 
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where  c  is  a  constant.  Now  we  see  from  (2.7)  that  all  of  the  functions  |u|,  |v|, 
|u»|,  jv'l  must  be  hounded  for  all  values  of  x  (since  p(x)  is  a  differenti- 
able  periodic  function  of  x)»  Therefore,  there  exists  an  upper  bound  for  Juv'-\ai*| 
which  is  independent  of  x.  According  to  (2.9)>  the  same  must  be  true  for  the  absolute 
value  of 

X 

I(x)  -  (     ^ 

However,  |l(x)|  -ooo  as  x  -^cd   since  u 

I(nn) 

Therefore  if  X  is  not  i-eal  we  cannot  have  a  solution  of  type  (2.7 ). 

Next  we  wish  to  show  that  there  exists  a  real  number  X  such  that  for  any 
X  <  X  the  solutions  of  (2.1)  are  unstable.  For  this  purpose,  select  a  X  such 
that  for  all  x 

X*  +  Q(x)  <  0. 

This  is  certainly  possible  since  Q(x),  being  periodic,  is  a  bounded  function  of  x. 
We  shall  show  that  if  X  <  X  then  y^(x,  X)  -ooo  as  x  ->a).  To  this  end  we  shaii 
write  (2.1)  in  the  form 

(2.10)  y-  -  D(x)y 

where  D(x)  -  -  X  -  Q(x)  >  0  for  all  x.  Since  y^(0)  -  1,  y^'  (0)  >  0  it  follows  from 

yj(0)  ■  0  that  y^(x)  >  0  for  all  srifficiently  small  positive  x.  Therefore,  if  the 

set  S  of  positive  zeros  of  y^(x)  is  not  empty  it  has  a  greatest  lower  bound  e  >  0. 

We  shall  show  that  e  does  not  exist  and  that  therefore  y ' (x)  >  0  for  all  x  >  0. 

For  this  purpose  we  observe  that  the  continuity  of  y,'(x)  implies  that  yHt)   ■  0. 

Now  we  have  from  (2.10)  that 

e 

(2.11) 


y{^(6)  -  2  j  D(x)  y^(x)y|(x)dx. 


Since  yj(0)  -  1  and  y|(x)  >Ofcr  0<x<ewe  have  y^(x)  >  0,   D(x)  >  0,   and 
y^(x)  >  0  for  0  <  X  <  e.     Therefore,   the  right  hand  side  of  (2.11)  is  positive 
whereas  we  had  assuried  that  yj^(6)  »  0.     This  shows  that  y'(x)  >  0  if  x  >  0.     There- 
fore yj^(x)  is  raonotonically  increasing  for  x  >  0  which  means  that  y,  (x)  >  1  for 
X  >  0.     Since  (2.10)  implies  that 


X 

y^^x)  -  2  I    D(t)yi(t)y{(t)dt, 
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we  see  that  yJ (x)  is  also  raonotonicaliy  increasing  with  x,   and  since  6  ■  yJ (x   )  >  0 

for  a  certain  x    >  0,  we  see  that 
o        ' 

y^(x)  >  1  +   (x-x^)5  (x  >  x^). 

Therefore  y,  (x)  ->oo  as  x  -»>od.  By  a  similar  argument,  we  can  show  that  yi(x)  >  1 
for  X  >  0»  We  have  thus  proved  incidentally 

Lemma  2«1»  If  X  <  X  ,  then,  f  cr  x  >  0, 
y^(x,X)  +  y^(x,X)  >  2. 

Since  this  result  implies  that  p,  ^   p_,  we  conclude  that  a  in  y. (x)  ■  A  e   p^ (x) 
+  B  e"  °^^p^(x)  cannot  be  real.  We  thus  find  that  if  X  is  complex  or  X  <  X  ,  then 
there  exists  no  solution  y(x)  of  (2.1)  which  is  of  type  (2.7)  with  real  a. 

We  shall  now  examine  closely  the  properties  of  the  functions  A  (X)  -  2  and 
^  (X)  +2,  where 

(2.12)  ^(X)  -  yj^(n,X)  +  y^(n,X). 

We  fist  note  that  ^  (X)  -  2  is  equivaloit  to  p.   ■  p^  ■  1  and  that  ^  (X)  «  -2 

is  equivalent  to  p^  ■  p^  ■  -1.     Hence,   if  ^  (X)  =  +  2,    (2.1)  will  have  a  solution 
of  type  (2.7)  with  real  a   (cf.  Floquet's  theorem).     Since  X  complex  or  X  <  X     implies 
that  (2»1)  has  no  such  solution  it  follows  that 

Lemma  2»2«     All  roots  of  the   equations  ^  (X)  -2  ■  0  and  ^ (X)  +  2  »  0  are 
real  and  >  X  , 

In  Section  2.2  we  shall  prove  that  both  ^  (X)  -  2  and  ^  (X)  +  2  are  entire 
analytic  functions  of  X  which  have  infinitely  many  zeros.     According  to  Lemma  2.2, 
all  of  these  zeros  are  real  and  greater  than  X  .     This  establishes  the  existence 
of  the  two  sequences   (2.2)  and  (2.3)  of  Theorem  2.1.     In  fact,  we  have  immediately 
from  Floquet's  theorem  that  the  following  assertion  is  true: 

Lemma  2»3.     Equation  (2.1)  has  a  periodic  solution  of  period  n  if  and  only  if 
A  (X)  •  2  and  a  periodic  solution  of  period  2n  if  and  only  if  ^  (X)  ■  -2. 

To  prove  Lenna  2.3»  we  merely  have  to  use  once  more  the  fact  that  the  condition 
p^   ■  p^  ■  1  of  Floquet's  theorem  is  identical  with  the  condition  /\  (X)  ■  2,   and  that 
p-   ■  Pp  ■  -1  is  equivalent  toy\  (X)  ■  -2. 


-  13  - 

Since  /\  (X)  is  an  entire  analytic  function,   the  limit  relations   (2.5)  are 
obviously  true.     We  now  turn  to  a  proof  of  the  ineqxwlities   (2.li).     For  this  purpose, 
we  need 

Lemma  2.U.     Let  n  be  a  root  of  the  equration  ^  (\)  -  2  ■  0  such  that  the 
derivative  ^  •  (X)  of  ^  (X)  with  respect  to  X  is  negative  or  zero  for  \  "  y.» 
Then  ^  "(X)  <  0  in  any  open  interval  m-  <  X  <  n^  in  which  ^  (X)  >  -2. 
Similarly,  let  n«  be  a  root  of  ^  (X)   +2-0  and  let  ^  '(ti')  >  0.     Then 
/^  '(X)  >  0  in  any  open  interval  ji«  <  X  <  m4  in  which  ^  (X)  <  2. 

Before  proving  Lemma  2.U,   we  may  observe  that  it  proves  both  the  inequalities 
(2.U)  and  the  assertion  that  the  open  intervals  (2.6)  are  intervals  of  stability. 
In  fact,  we  see  from  Lemma  2.1  that  if  X  <  X  ,  then  ^  (X)  >  2.     Among  the  infinitely 
many  real  zeros  of  the  function  A  (X)  -  2  there  must  be  a  smallest  one  which  we 
call  X  .     We  shall  prove  later   (see  Lesamia  2.6)  that  A    '(X  )  >  0.     Therefore  Lemma 

O  *""*    o 

2»l4  shows  that,  for  X  >  X  ,  A  (X)  must  be  a  decreasing  function  until  ^  (X)  ■  -2. 
This  must  actually  happen  for  a  certain  X  ■  X»  >  X  ,  since  ^  (X)  +  2  has  infinitely 
many  real  zeros  without  a  finite  limit  point.  Now, either  ^  '(XJ)  •  0,  or 
^  '(X")  <  0.  If  /2i  '('^^  "  0»  X»  is  a  double  root  of  ^  (X)  ♦  2  »  0  and  will  be 
listed  as  X*  and  as  X».  (According  to  Lanma  2.5  proved  below  ^  (X)  +  2  ■  0  cannot 
have  roots  of  multiplicity  higher  than  two.)  ^  '(XJ)  -  0  implies  (of.  Lanma  2.U) 
that  ^  (X)  increases  for  X  >  X'  ■  X»  until  it  reaches  the  value  2.  On  the  other 
hand,  if  ^  '(X*)  <  0,  then  ^  (X)  <  -2  for  X»  <  X  <  X »  ,  where  X"  is  the  smallest 
zero  of  ^  (X)  ♦  2  which  is  >  X«.  Since  ^  (X)  <  -2  in  the  interval  (X«,  X')  this 
is  an  interval  of  instability  for  the  solutions  of  (2.1)  (cf .  the  stability 
test  of  Section  1.2).  Now,  ^  '(XJ)  >  0  as  can  be  seen  from  the  fact  that  ^  (X)  < 
y\  (X»)  for  all  X  <  X«  and  sufficiently  close  to  \X»     Using  Lecma  2.U  we  may  there- 
fore conclude  that  ^  (X)  is  an  increasing  function  of  X  in  any  interval  X '  <  X  <  X 
in  which  ^ (X)  <  2.  The  largest  interval  of  this  kind  is  the  interval  (X«,  X^), 
where  X^  denotes  the  smallest  root  of  ^  (X)  -  2  which  is  >  X«.  The  stability  test 
of  Section  1.2  shows  that  this  interval  is  an  interval  of  stability  for  the  solutions 
of  (2.1). 

Continuing  in  this  manner  we  find  that  the  inequalities  (2.1;)  hold  and  that 
the  open  intervals  (2.6)  are  the  only  intervals  of  stability  for  the  solutions  of 
(2.1). 

We  now  prove  Lemma  2.14.  For  this  purpose,  we  introduce  the  following  notations: 

§^  yi(x,X)  -  z^(x,X),    1^  y2(x,X)  -  Z2(x,X), 
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1^  y|(x,X)  -  z^(x,X),         1^  y^(x,X)  -  z^(x,X) 

where  obviously  zJ  ■  ^  z.   and  ^n  '  "Sx  ^2*     ^^°»  *"  shall  write   y^,    Jj?^,  )f  ^,  i^^ 
respectively  for  y,  (n,X),  ypC^jX),  y^(n,X),  y'(n,X),     As  before,  we  shall  write 
/\for  "yj-,   *y)X  s^nd  ^'  for  the  derivative  of  ^  with  respect  to  X,   that  is: 

^  -  Zj^(ji,X)  +  z^(n,X). 

Next,  we  shall  derive  the  formula 

n  n  n 

(2.13)   ^'(X)  -  (>^i  -  ^p  [  y3^(x)y2(x)dx  -  i^  g  (  ^1^^^^  *  ?  1  f  A^^"^^^' 

'o  'o  *0 

To  prove  (2«13),  we  differentiate  (2.1)  with  respect  to  X  and  we  obtain  (for 
y  «  y  and  y  "  yp  respectively) 


(2.1U) 


z{'+  (X  +  Q)z^  -  -y^ 


z^'  +  (X  +  Q)z2  «  -yg 


The  general  formula  for  the  solution  cf  an  inhomogeneous  linear  differential 
equation  of  the  second  order  in  terms  of  the  solution  of  the  homogeneous  equation 
yields 


5^(x)  -  y^(x)  j  y2(t)y^(t)dt  -  y2(x)  (     y^(t)dt 


o  o 

X  x 


(2.15) 


z{(x)  -  y{(x)  I  y2(t)y3^(t)dt  -  y^(x)  j  y^(t)dt 
/o  'o 

X  X 

Z2(x)  -  y^(x)  [  y2(t)dt  -  y2(x)  [  y^(t)y2(t)dt 


(    y2(t)dt  -  y2(x)  f 
'  0  'c 

X  X 

2^(x)  -  y W  y2(t)dt  -  y^(x)  f  yj^(t)y2(t)dt  . 


The  functions  z^   and  z,  in  (2.15)  are  those  solutions  of  (2.1U)  which  satirfy  the 
initial  conditions  z^(0)  -  z|(0)  ■  0,  Z2(0)  -  z^(0)  »  0.  Since  z^  -  (^)yj  etc., 
and  since  the  initial  conditions  for  y^  and  y-  are  independent  of  X,  tie  solutions 
(2.15)  are  the  coirect  ones.  Equation  (2.13)  follows  immediately  from  (2.15)  by 
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putting  X  ■  n» 

Since  the  Wronskian  y-,y^  -  y^y-!  ■  1  for  all  x,   we  find  for  x  ■  n  that 

(2.36)  7l72-   nz   Vl'^- 

Hence 

Putting  sgn   v^  »  ■  +  1  if    >^{  >  0,   sgn   wJ^-lif    V7{<0  and  sgn   y?  '   -  0  if    >^»  ■  0, 
and  assuming  that   x' •   ji^  0,  we  find  from  (2,13): 


^•(\)  -  sgn  y^^^   fh/r^{|    Y2  *  ^^^\{ 


ni-nz 


(2.17) 


Equsticn  (2.17)  shows  that  ^'(X)  has  the  same  sign  as  )^  •  in  any  interval  in  which 
y«J  /  0  and  ^  <  U.  Consider  now  a  value  \i.  of  X  such  that  A   (p.)  ■  2  and 
^•(n.)  <  0.  We  wish  to  establish  the  fact  that,  for  a  sufficiently  small  5,  ^  (X) 
is  decreasing  in  the  interval  tx  <  X  <  X  ♦  6.  If  ^'(p.)  >  0,  this  is  obvious.  Assume 
now  that  ^  {\i.)   -  2,  £^{\i)   ■  0.  In  this  case,  according  to  (2.17),  we  must  have 
yi  1  (p.)  -  0.  Since  we  also  have 

A^  -  i^  •  ( yi  -  >?i)^  *  ^  ^{  h  "  ^' 

we  find  that   >?  A\i)  -    y7Xiii)  -  0  and,    sxnce  *^  ■i'^  i  I2  ^i  '  ^*  ^^  '^^^® 

v^j^(ti)  -  >^A(m-)  ■  !•  Then (2.13)  reduces  to 

n 

'o 

and  therefore  ^' (p)  ■  0  implies  ^^(M')  "  0. 
Now  we  shall  compute 

for  X  »  ti,  where  [l  is  such  that 
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(2.18)  y^[i^)  -  vi^iv)  -  0,       y^(ii)  -  yi'^iii)  -  1. 

We  shall  do  this  by  differentiating   (2.13)  with  respect  to  X  and  by  using   (2»15) 
for  X  ■  n  in  order  to  obtain  the  derivatives  with  respect  to  X  of    >?,(X),   etc.,  for 
X  ■  ii»     A  straightforward  computation  shows  that,   if   (2»l8)  holds, 

/n  \  n  n 

I       yT(x)y^(x)dx  J^  -  2j     y?(x)dx  (    y^(x)dx     . 


Since  y, (x)  and  y^Cx)  are  linearly  independent  functions,  we  see  from  an  application 
the  Schwartz  inequality  to  (2.19)  that 


II 


(2.20)  ^   (^)  <  0  . 

Thus  A  (X)  is  again  found  to  be  decreasing  in  an  interval  |i.  <  X  <  |i-»-5.  Assume 
now  that  2»k   is  false*  Then  there  would  exist  a  smallest  number  li  >  ji  such  that 
^  (X)  <  0  for  M-  <  ^  <  M-  but  ^  (ji  )  -  0  although  ^(n  )  >  -2.  We  would  then 
have 

(2.21)  A^/)  -  U  -  (  >^1-  y^2^^  *k   >^i  ^2  *  ° 

and  therefore  Y)^  %  "^  ^  ^^"^  X  ■  p.  ,  But  then,  )9  ••  (m^  )  ¥  0  and  so,  according  to 
(2.17),  /\(m.  )  /  0, which  produces  a  contradiction.  This  proves  Lemma  2.U  in  the 
case  where  /\  (X)  ■  2. 

If  ^  (X)  «  -2,   the  proof  is  almost  literally  the  same.     Incidentally,   our  proof 
of  Lenima  2»U  shows  that  the  following  is  truei 

Lemma  2.S.     The  roots  of  the  equation 

^M  -  U  -  0 

are  either  simple  or  double  roots.  If,  for  a  particular  value  of  X  ■  ^l, 

(2.22)       ^i\x)   -  U,   ^'(n)  -  0, 

then^  (n)  <  0  if  ^(m^)  ■  2  and  ^   (n)  >  0  if  ^-  -2.  Necessary 

and  sufficient  conditions  for  ^  (p.)  -  U  and  ^  (^)  to  vanish  simultaneously 

are 


(2.23)     y^^M  -  yi^M  -  yi[i\i)  •  yi^iik)  -  o, 
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In  order  to  complete  the  proof  of  Theorem  2»1  we  need 

Lcpgna  2>6»  Let  X  be  the  smallest  root  of  the  equation  /\  (X)  -  U  ■  0. 
Then  X  is  a  simple  root  and  ^  (X  )  <  0« 

Proof*  We  know  from  Lemma  2*1  that  /\(X)  >  2  for  X  <  X  •  Therefore  \  '  \    cannot 

be  a  maxiHum  of  A(X)  if  ^(X  )  ■  2.  But  Leinna  2.5  shows  that  /\(X)  would  have  a 

maximum  at  X  ■  X  if  A  (X  )  ■  0»  This  is  a  contradiction,  and  therefore  Lemtia  2.6 
o    ^--*   o 

is  true* 

A  comparison  of  Lennna  2«5  and  of  Floquet's  theorem  in  Chapter  I  shows  immediately 
that  we  can  supplement  Theoron  2»1  by  the  following 

Corollary  2.1»     Hill's  equation  (2»1)  has  two  linearly  independent  periodic  • 

2 
solutions  of  period  n  or  2n  if  and  only  if  the  equation  ^  (X)  -  U  ■  0  has 

a  double  root. 

The  stability  test  at  the  end  of  Section  1.2,   Corollaiy  2.1  and  Lemma  2.5  show 
that  the  solutions  of  (2«1)  are  stable  for  ^  '  \  ^i   (or  X  ■  X"   ^. )  if  and  only  if 
\J^^.1  (or  ^n+1^  ^^  *  double  root  of  ^  (X)  -  U,  that  is,  if  and  only  if  X.  ^,  • 
>2n+2   ^°^  ^n+1  '  Mn+2^*     "^^^  proves  Theoran  2.1  completely. 

2.2.     Analytic  properties  of  the  discriminant 

The  function  ^(X)  as  defined  ly   (2.12)  will  be  called  the  discriminant  of 
Hill»s  equation  (2.1).     In  this  section,  we  shall  prove  a  result  which  we  used 
already  in  Section  2.1  and  which  we  may  state  as  follows: 

Theorem  2.2.     The  function 

^M  -  yi(n,   X)   ♦  y2(n,X) 

is  an  entire  analytic  function  of  the  complex  variable  X.     Its  order  of 
growth  for  |x|  ->  oo  is  exactly  ^  }  i.e.,  there  exists  a  positive  constant 
M  such  that 


(2.2U)  |/\  (X)|    exp(-M  /-fXf  ) 

is  bounded  for  all  X  and  a  positive  constant  m  such  that,  X  real 
and  X  — >  -  00  implies 

(2.25)       \/:^W\    exp(-m  /"PT  )  ->  oo  . 
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Corollary  2»2.     The  functions  ^(X)  ♦  2  and  /^(X)  -  2  have  infinitely 
many  zeros. 

Note  that  Corollary  2*2  follows  immediately  from  Theorem  2*2  which  permits  us 
to  conclude  that  ^(X)  +  2  and  ^(X)  -  2  are  functions  of  order  of  growth  r-  . 
According  to  a  well  known  theoran  on  entire  functions,    (see  Nevanlinna,   1936  or 
Titchraarsh,  1938)  any  entire  function  of  order  of  growth  «•  has  infinitely  many  zeros. 

In  order  to  prove  Theorem  2.2,  we  apply  Picard's  method  of  iteration  to  the 
differential  equation  (2.1).     Let  co  "  /T  and  let 

sinojx 

U      ■   cos  ftBC,  v^   •   

O  '  0  0) 

and  define  u  ,   v    recursively  for  n  ■  1,   2,   •••by 

X 

(2o26)  u^(x,co)  -  -  I  J     sin  a)(x-4)  Q(4)  Vl^^^*^ 

o 

X 

(2.27)  v^(x,(»)  "-11     sin  ®(3^«)  Q(4)  Vl^^^'^     * 

'o 

Then  ^ 

(2.28)  yi(x,©  )  -  5~  u-(x>a)) 

n"o 

(2.29)  y,(x,co'')  -^v^(x,co) 


n"o 


and 


(2.30)  ^(X)  -^(co^)  -^^I  PnC"*'^)  +  v^(«>«)] 

n*o   l-  — 1 


where  v»(3C,aj)  ■  cos  cox  and 


(2.31)                v^(x,co)  -  - 

cos  «(x-4)  Q(«)Vj^ 
o 

It  is  easy  to  see  that  for  all  real  values  of  x  >  0: 

K  i*"' 

.  \\  <-'''=' 

n-1' 
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Mow  let  M  be  a  positive  constant  such  that,  for  all  real  values  of  x,  |Q(x)|  <  M  • 
Then  we  find  by  induction  from  (2.26),  (2.27)  and  from 

|sin  (o(x-4)|  <  |«|  (x-4)  el'^l  ^""^^  (0  <  4  <  x) 
that 

(2.32)  |u^(x,a.)|  <  el"l  V  x^)"  /(2n): 

(2.33)  |v^U,co)|  <  X  e'^^I^M*  x^)^  ^  (2n*l):  . 
Equations  (2.33)  and  (2.31)  show  that 

(2.3U)  |v;(x,a))|  <  •'"^'^(M*  x^)"  y/  (2n): 

and  we  therefore  conclude  from  (2.30)  that 

(2.35)  lA(co^)|  1  2  el'^l"  cosh  {^   n)  . 

This  proves  that  the  expression  (2.2U)  is  bounded  if  we  choose  M  ■  n.  In  order 

to  complete  the  proof  of  Theorem  2.2,  we  have  to  prove  (2.25).  For  this  purpose, 

* 
we  may  assume  that  Q(x)  <  -1  for  all  x.  Othei*wise,  we  could  replace  Q  by  Q-M  -1 

and  X  by  X  +  M  +1  without  changing  the  differential  equation  for  y.  Putting 

^  ■  i(?  -  CO,  where  (J  is  real  ar^l  positive  and  j?-«»ao  asX.^-QD,we  find 


1  0x   1  ax/2 
%  >  2  ^   ^  ?  * 


/-4 
sin  co(x-€)  ,  sinh jr(x>»^)  _  1  f   g^*  ds  > 
(0         i?        2  /         — 


'(x-?)/2 


From  the  preceding  inequalities  and  from  (2.26)  we  find  by  induction  with  respect 
to  n  that 
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(2.36)  u^(x,icr)  >^  6^^/2(^/2)2'^    /   (2n):   . 

By  an  even  simpler  argvment,  we  can  prove  that  v' (x,i(2f)  >  0,   md  therefore  we  find 
fom   (2.36)  and  (2.30)  that 

(2.37)  ^M  -  ^i-^)  >  ZI  %(«»ii2l)  >  \  e^"/2  ^^^^  (^/g)     . 


This  proves  that  (2.25)  is  true  for  any  »  between  0  and  b/2. 

Obviously,  all  the  functions  u  and  v'  are  entire  analytic  functions  and  X 
and,  since  we  have  shown  that  their  sums  converge  uniformly  in  any  finite  part  of 
the  X-plane,  it  follows  that  ^(X)  -  yj(n,X)  ♦  y2(n,X)  is  also  an  entire  function 
of  X.  The  inequalities  (2.35)  and  (2.37)  establish  the  truth  of  our  assertion 
about  the  order  of  growth  of  ^(X)  and  this  completes  the  proof  of  Theorem  2.2. 

For  later  purposes  we  note  here  a  result  which  can  be  proved  by  exactly  the 
same  method  as  Theorem  2.2. 

Theorem  2.3.  Let  7(x,X)  be  any  real  solution  of  Hill's  equation  (2.1)  with 
initial  conditions  independent  of  X.  Let  x  be  fixed  and  real,  and  consider 

y(x,X)  and  y'(x,X)  as  functions  of  X.  Then  the  order  of  growth  of  these 

J——  -  —^ 

two  fiinctions  of  X  is  at  most  «•  • 


The  following  theorem  gives  an  idea  about  the  asymptotic  behavior  of  ^(X)  for 
positive  values  of  X. 

Theorem  2.U.  The  absolute  value  of  the  function 


A     r^C^^)  -  2  cos  n^T  I 


is  bounded  for  all  real,  positive  values  of  X. 

Proof.  If  X  -  oo^  >  0,  we  find  from  (2.26)  and  (2.27)  (by  induction)  that  for 
X  >  0  and  n  ■  1,  2,  3,  •••« 

|u„(x,a>)|  <  (M*)"xV^/n: 

|vJx,co)|  <  (MVxV^"^/n: 
The  estimate  for  |v  |  together  with  (2.31)  shows  that 

Iv;^(x,co|  <  (H*)"  x^co"" /nj 
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Therefore  we  see  from   (2.28)  and   (2,29),   that,  for  x  >  0, 

(2.38)  |y^(x,co^)  +  y2(x,co^)  -  2  cos  (ox|  <  2  exp(7M*  /«)  -  2 

and,   in  particular,  for  x  ■  n: 

|^(X)  -  2  cos  lA  n|  <  2  exp  ~  -  2  . 

Since  for  any  real  t  >  0 

t 

exp  t  -  1  -    exp  3  ds  <  t  exp  t 


■( 


we  find  that 


It  It 

(2.39)  |A(X)  -  2  cos  «  J/T  I  <  22L  exp  :2L 

which  proves  Theorem  2*U. 

Obviously,  Theorem  2«U  is  true  for  all  bounded,  square  integreble  functions 
Q(x)«  In  order  to  obtain  better  results  pertsdning  to  the  asymptotic  behavior  of 
/\(X)  for  X  >  0,  X  ->  00,  we  shall  assume  that  Q(x)  can  be  expanded  in  a  Fourier 
series 

(2.U0)  Q(x)-2I'   gn*^^°^' 

n  -  -  00 

where  the  prime  at  the  staoiation  symbol  indicates  that  the  sum  is  extended  over 

values  of  n  y  0  only.     This  does  not  impose  restriction  on  Q(x),   since  a  constant 

term  in  (2.1j0)  could  be  combined  with  the  parameter  X  in  Hill's  equation  (2.1). 

Equivalently,   we   could  say  that  Q(x)  has  been  normalized  so  that 

n 

(2.1il)  f     Q(x)  dx  -  0 

0 

Since  Q(x)  is  real,   the  constants  g    in  (2.1iO)  satisfy  the  condition 
(2.U1*)  g.^  .  g„ 

for  all  n,  ^ere  a  bar  denotes  the  conjiogate  complex  quantity.     We  shall  always 
assume  that 

+  00 

(2.U2)  2Z      l8j  >~ 

n»-oo 
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and  In  many  cases  that 


(2.U3)  lim     "^  8n  "  0| 

n  — >  +00 


(this  will  be  true  if  Q(x)  has  a  continuous  second  derivative). 
If  we  put 

(2.UU)  /:^(X)  -  u^(n,  /T)  *  v^(n.  A)  , 

where  u  ,  ▼  are  defined  by  (2»26)  and  (2.27),  then 
n'  n 

Obviously,  /\_(X)  is  a  homogeneous  form  of  degree  n  in  the  infinitely  many  variables 


g  ,   that  is 

(2.145) 

/^M  - 

r 

h""*  ^n""" 

c(/p...,  /^)  g/  ...g^ 


where  we  may  sxm  without  restrictions  if  we  assume  that  g  "0.  Now  we  shall  prove: 
Theorem  2.5«  Let  oo  ■  yT  .  Then 

c(£i,...,i„)  -  A(a>)  cos  no)*  B(a))  5iDJH  , 


where  A{<a)  and  B(<a)  are  even  rational  functions  of  co  such  that 

(1)  In  each  of  the  functions  A(co)  and  B((o)  the  degree  of  the  denominator 
exceeds  the  degree  of  the  numerator  by  at  least  n. 

(2)  The  poles  of  A(co)  and  B(a))  are  at  some  of  the  points  cd  -  0  and 

where  1  <  r  <  «  <  n. 

Proof  of  theorem  2.$.  We  shall  use  the  expression  for  ^   in  terms  of  u  (n,\) 

and  v^(n,X),  However,  the  expressions  of  u  and  v»  as  given  by  (2.26)  and  (2.31) 
involve  n  integrations,  and  direct  evaluation  of  these  integrals  leads  to  results 
which  are  cumbersome  in  form.  We  shall  therefore  make  use  of  the  fact  that  the 
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integrals  in  (2.26)  and  (2.31)  are  convolution  integrals  the  Laplace  transforms 
of  which  can  be  easily  computed.  We  find  for  the  coefficient  of 


g;  gf      ...  g/ 


in 


j   e"P^  [u^(x,X)  ♦  vi(x,xjj  dx 


the  expression  n 

2p.  2iTt 
^  v«i 

(2.1i6)    k(j^^,...,/^,p) 


p^p^  n    p  *   (p-2i/^  -  ...  -  2ii^)^ 


By  an  application  of  the  inversion  f onoula  for  the  Laplace  transformation  we  find 
that  c(f^,..»/  )  is  the  sum  of  the  residues  of 

U  -  9"P  k(j^^,...,4,p). 

Obviously,  the  poles  of  U  are  located  at 

p  -  ico,  ico  +2i/p  i«  +  2i/^  +  2i/2»... 

and  Theorem  2.5  follows  now  from  an  inspection  of  U  in  the  neighborhood  of  its 

poles  and  from  the  remark  that  c(/^, ..,j^   )  must  be  an  even  function  of  co  (since 

it  is  an  entire  function  of  X)« 

Since  the  product  g£  ...g^    does  not  change  if  we  permute  the  variables 

1  n 

g/  f*tSf  i  we  find  that     /\     (X)  can  also  be  written  in  the  form 
1  *n  " 

Where  -oo  <  / 1  1  i^p  •**  -  ^n  "^  ^'     "^^^  coefficients  g  in  (2.U7)  are  defined  by 
the  relation 

(2.U8)  ta^,    ...,/„)   -^C(K^,    ...,    K^), 

Where  the  sum  in  (2.U8)  is  to  be  extended  over  all  distinct  sets  of  integers 
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K,,  ••»,  K  which  can  be  made  to  coincide  nith  the  set  /^,  ...,  /^  by  a  suitable 
permutation  of  the  indices  1,  ••»,  n.  Concerning  the  coefficients  y  we  have 

Theorem  2*6.  The  coafficients  y  defined  by  (2«U8)  satisfy  the  relations 
whenever 

H^  *  ...  +  /„  /  0. 

Proof*     If  we  replace  Q(x)  in  (2.1)  by  Q(x  +  5^),  where  flf  is  a  constant,   the 
resulting  differential  equation  will  have  the  same  discriminant  ^(X)  as  the 
original  one.     In  fact,  it  is  clear  that  if  Q(x)  is  replaced  by  Q(x  +  0),  Hill's 
equation  continues  to  have  periodic  solutions  for  the  same  values  of  X  and  in 
each  case  their  number  is  unaffected  by  the  change.     Since  ^(X)  -  2  is  an 
analytic  function  of  X  whose  order  of  growth  is  j-,   ^(X)  -  2  is  determined  by  its 
zeros  up  to  a  multiplicative  constait.     Assume  now  that  ^(X)  -  2  belongs  to  the 
function  Q(x)  and  that  ^*(X)  -  2  belongs  to  the  function  Q(x  ♦  0).     Then 


^(X)-2  -  c[^*(X)-2]  , 


Where  C  is  a  constant.  On  the  other  hand,  ^(X)+2  and  ^  (X)+2  also  have  the 
same  aeros  (with  the  same  multiplicity),  and  therefore 

^(X)+2  •  c'j^*(X)*2l. 

By  subtracting  the  first  of  these  aquations  from  the  second  we  find  that 

(c'  -  C)  ^*(X)  -  2C'  +  2C  -  1*. 

Therefore  A  (X)  would  be  a  constant  unless  C  ■  C  ■  1.  However,  Theorem  2.U 
shows  that  A  (X)  can  never  be  a  constant.  Hence  ^(X)  ■  ^  (X)  for  all  X. 

It  is  obvious  that  /^*(X)  arises  from  ^(X)  by  substituting  g^e  ^^   for  g^, 
and  therefore  the  coefficient  of  g/   ...  g^  in  ^  will  be 

g(^l,  ...,  i„)  e^il^, 

where  L  ■£,  +  ...  *  I    .  Assume  now  that  all  of  the  g^  are  zero  except  for 
g/  »  Zt  *    •••»  g;  •  The  sane  argument  which  we  used  in  showing  that  ^(X)  is 
an  analytic  function  of  X  (Theorem  2.2)  can  be  used  to  show  that  ^  depends 
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analytically  on  g/i  ^   ...,    %o    and  the  same  is  true  for  ^  .     Hence  ^(>^)  -  /\  (X), 

implies  that  the  coefficients  of  g/!'     .•»  gy    in  both  functions  must  be  the  same, 
so  that  for  all  real  values  of  <$x 

If  L  ^  0,   this  is  possible  only  if  y  "  0»  as  stated  in  Theorem  2«6.     We  may  note 
here 

Corollary  2.6>     The  first  terms  in  the  expansion  (2.1tU)  are 
^^  (X)  -  2  cos  n  A,  Zli  ^^^  "  ^ 


.      /^x       71  sm  n  /A.      ^— •       '"n' 


2  A 

Of  these  relations,  the  first  one  is  trivial,   the  second  one  follows  fran  Theorem 
2e6,   and  the  third  one  is  a  result  of  a  somewhat  tedious  computation. 

2.3     Infinite  Determinants 

Hill,   1886,   used  infinite  determinants  for  the  investigation  of  the  character- 
istic values  of  X  in  (2.1).     Whit taker  and  Watson,  1927,  showed  that  the  value  of  Hill's 
determinant  can  be  expressed  in  terras  of  a(X).     In  this  section  we  shall  reproduce 
the  results  of  Hill  and  Watson  and  supplement  them  by  some  relations  of  the  type 
obtained  by  Whittaker  and  Watson. 

We  shall  write  a  determinant  in  the  f oiw 

'       k 

where  n  and  m  vary  over  all  integers  from  k  to/.  In  particular,  we  shall  consider 
the  determinants  where  k  •  -oo,  i  •  oo  or  where  k  "  0,  j^«  oo.  These  we  shall  call  two 
sided  infinite  and  one  sided  infinite  determinants,  respectively. 

We  shall  always  use  the  first  subscript  n  to  denote  the  rows  and  the  second 
subsciript  m  to  denote  the  columns  of  the  determinant. 

We  shall  say  that  the  infinite  determinants 

II  a   r  ,  ||a   II  °° 
"  n,m"  o  *  "  n,m"  -oo 

exist  or  converge  if  the  limits 
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lim       II  a       11^      ,     Urn       ]|a       1(^ 

exist.     The  value  of  the  limit  is  then  called  the  value  of  the  determinant.     We 
shall  not  be  concerned  here  with  a  general  theory  of  infinite  determinants;  instead, 
we  shall  introduce  a  special  class  of  infinite  determinants  arid  indicate  briefly 
the  proof  of  a  few  theorems  which  hold  for  this  class. 

We  shall  say  that  a  determinant  is  of  Hill's  type  if  it  satisfies  the  condition 

*>—     '   n.Bi        n.m 
n,m       ^  ' 

where  6         "1  for  n  ■  m  and  5         "0  otherwise,    and  where  the  sum  in  (2.U9)  is 
n,m  n,m  * 

to  be  taken  over  all  values  of  n  and  m*     Obviously,   every  finite  determinant  is 
of  Hill's  type.     We  shall  show  now: 

Theorem  2.7.     An  infinite  determinant  of  Hill's  type  converges. 

It   suffices  to  prove  Theorem  2.7  in  the  case  of  a  determinant   lla       11^  • 

"^n,m"  o 

According  to  a  theorem  due  to  Hadamard  (see  Hardy,  Littlewood,  Polya,  193U)  the 
absolute  value  of  the  square  of  any  finite  determinant 

ll^n,mllt 
does  not  exceed  the  value  of  the  product 


V 

ii"k 


(i  ".-'^J 


From  this  we  derive  the  following 

Lemma  2.7.     Let  ||  a       11        be  a  determinant  of  Hill's  type.     Let 


a'       -a        if  n  j'  m 

n,m         n,m 


a»        -  a        if    a  >  1 

n,n  n,n         '   n,n'   — 


a«       -  1  if     a  <  1, 

n,n  '   n,n'         ' 


and  let 
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03  /OO 

H  -<     IT        Z     Ia' 
n-0     V"*^ 


Then  H  <  oo,   and  the  absolute  value  of   any  finite  subdetemiinant  of 
II  a^  ^11      is  at  most  equal  to  H. 

Proof.   The  only  difficulty  is  to  show  that  H  <oo»     For  this  purpose,   let 

e     =T'    |a       I 

m  * 

where  the  sum  is  taken  over  all  m  /  n.  Also,  let 


Y  "  a'    -  1 
'n   '  n,n' 

It  follows  from  (2.1i9)  that 

00  00 


n-C  "       rK)  " 


Finally,  let 


00 


1 1/2 


^n-^^KJ  - 


m*0 


Then 


V    <  |a'   1^  -  1  ♦  f  r'  la   II 
"-  '  n,n'        Uj.  I*n,m"y 


|a'  ?■   -  1  ♦  c^<  A|a   -  1|  +  6^ 
n,n'        n  —  '  nn    '    n 


where  A  ■  max  a    +  1, 
'  nn' 


In  view  of  (2.U9) 


CO 


5^  a   -  1  <  00  . 

^—x   '  nn    ' 
n^ 


Also,  we  know  that  the  convergence  of  a  sum  of  positive  terms  implies  the  convergence 
of  the  sum  of  the  squares  of  the  terms*  Hence 

00  r. 


S""' 


00 


and 


H 


fr  (1  ♦  P^) 

n"0 


-1 1/2 


<    00 


This  proves  Lemma  2.7» 

Now  we  can  prove  Theorem  2.7.     We  can  expand 


%mlir       •  H^Y 
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in  terms  of  the  elements  of  the  last  row  and  theii*     3ubdeterminants     which  can 
be  majorized  by  the  quantity  H  defiijed  in  Lemma  2.7.     The  result  is 

where   \   %\  <  ^Ui*     Since   |D^|   <  H,  it  follows  that 


00 

z 

m-0 


«^       l^/+l,m  -  SM,ml 


Therefore 

and 

exists.     This  proves  Theorem  2«7» 

The  last  result  about  infinite  determinants  which  we  need  is 

Theorem  2.8.     Let   11  a       ll    be  an  infinite  determinant  of  Hill's  type, 
—  i    -       — —   "    n^m'  ■  ■»  ' 

and  assume  that  there  exist  numbers  x,  not  all  of  which  vanish  such  that 

I  m  — 

|x  I  <  M  (M  fixed)  for  all  m  and 

^-    n.i 
m   ' 

for  all  n.  Then 


X  -  0 
im  m 
m 


II  a   II  -  0. 
' '  n,m ' ' 

Proof.  Since  the  set  of  subdeterminants  of  ||  a   ||  is  bounded,  it  follows  that 
for  each  n 

II  a   ||x  -  0. 
' '  n,m "  n 

This  relation  is  obtained  in  exactly  the  same  manner  in  which  the  corresponding 
relation  is  obtained  in  the  case  of  a  finite  system  of  linear  equations.  The 
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inequality  (2.U9)  and  the  condition  |x  |  <  M  guarantee  the  absolute  convergence 
of  all  sums  involved.     This  proves  Theorem  2.8. 

We  shall  now  express  the  discriminant  /\(X)  of  Hill's  equation  in  terms  of 
an  infinite  determinant.     For  this  purpose  we  shall  write   (2,1)  in  the  form 

(2.50)  y  +  (I:   g y^"^ )  y  -  0 

\n--oo  "    / 

where  X  ■  g  and  Q(x'*  is  givam  by  (2.U0).  We  know  from  Floquet's  Theorem  that 
(2.50)  has  a  solution  /  0  of  the  type 


(2.51) 


iax  /  V 
7  •  e    p(x) 


where  p(x)  is  a  function  of  period  n  and  where 

2  cos  na  •  y^(Ti)  ♦  ygCn). 

If  Q(x)  is  sufficiently  smooth,  e.g.  if  ^  |g  |  <  cd  ,  the  function  y  in  (2.51) 
can  be  expanded  in  a  twice  termwise  differentiable  series 


CD 


i(a+2n)x 


(2.52)  y(x)  -  Y.      Pn' 

n=-aD 

and  the  left  hand  side  of   (2.50)  takes  the  form 
(2.53) 


f    C  ei(^*2n)x 
^—   n 


-00 


Since  (2.53)  must  vanish  identically,  we  have  C  •  0  for  all  n.  If  we  write  C 
explicitly  in  terms  of  the  p  and  g  ,  we  have,  for  -oo  <  n  <  oo. 


(2.5U)  I:   fg^^  -  (a*2n)2  6  1  P„  -  0, 

or,  after  multiplication  by  g  -  (a+2n)  ~  : 


(2.55) 


00 

z 

m=  -  oc 


\-i 


m 


X-(a+2n) 


2     n,m 


m 


where  we  have  replaced  g  by  X  and  where  now.  Just  as  in  (2.U1  ) 


(2.56) 


^n-m  ■  Vn  '   ^^  "  °  ' 


Obviously,  the  determinant 
(2.^7)  D(c,X) 
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'n-ffi 


T  -^  6 

X-(a+2n)      '  I 


oo 


00 


converges  if  ^  |g  |  <  oo,  except  for  such  values  of  X  and  a  for  which  one  of 

2 
the  denominators  X-(a+2n)  vanishes. 

It  is  easy  to  see  that 

Lemma  2»8«  D(a,X)  regarded  as  a  function  of  o  is  single  valued  and 
analytic  for  all  values  of  a  other  than  the  values 


(2.58) 


■  +  A  -2n,   n  ■  0,  +1,  +2,  see  , 


at  which  the  function  may  have  poles.  If  X  j^  0,  these  poles  are  (at 
most)  of  order  one.  D(a,X)  is  periodic  (in  a)  with  period  2  and  for 
a  -i-   ioo,  D(a,X)  ->  1« 

The  proof  of  Lemma  2.8  is  mostly  routine  except  for  the  statement  about  the 
periodicity  of  D(o,X).  This  follows  from  the  remark  that  D  remains  unchanged  if 
we  replace  a  by  a+2  and  at  the  same  time  replace  n  by  n-1  and  m  by  m-1.  Since  both 
n  and  m  run  from  -co  to  +00,  the  same  is  true  for  n-1  and  m-1,  and  therefore  D  does 
not  change  if  we  replace  a  by  a+2.  The  rest  of  the  proof  of  Lemma  2.8  is  left  to 
the  reader. 

Since  the  residues  of 


n-m 


X-(a+2n)' 


at  the  values  a  -  yT  -  2n  and  a  ■  -  /T  -2n  add  up  to  zero,   it  follows  from  the 
periodicity  of  D(a,X)  that   (for  \^0)  all  of  its  residues  have  the  same  value  K 
for  a  «  jA  -  2n  (independent  of  n)  and  the   value  -K  for  a  ■  -  yT  -  2n. 
Therefore 


(2.59) 


E(a)  -  D(a,X)  -  K  • 


ctg  n.(a  -  A)  -  ctg^a+A  )' 


is  an  entire  function  of  a  with  period  2.     Since  E(a)  is  bounded  in  the  strip 

-]<Rea<l     , 


E(a)  is  a  constant  E.  We  wish  to  determine  E  and  K.  By  letting  a  ->  ioo,  we 
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find  from  Lemma  2.8  that  E  -  1.  We  cannot  determine  K  explicitly,  but  we  can 
express  it  in  terms  of  D(0,X)  by  putting  a  ■  0.  The  result  is 

(2.60)  K  -  I  |l  -  D(0,X)j  tg  J  /I  . 

We  can  now  apply  Theorem  2.8  to  (2.55).  The  infinitely  many  equations 
(2.55)  are  not  always  the  equivalent  of  (2.5U),  since  a  may  have  one  of  the 
exceptional  values  (2.58).  However,  we  may  multiply  (2.55)  by 

(1  -  LzJL.)  (1  ♦  °  -^  >^  )  -  (2n-^a)  -  X 
2n  2n         Un 

for  n  /  0  and  by  a     -  X  for  n  ■  0.     The  determinant  of  the   resulting  system  must 

vanish  if  not  all  of  the  p     vanish.     Since 

m 

"   (a^  -  X)  TT     (1  -  ^^)    (1  *  ^)   -  sin  i(a-A)  sin  }{a./T  )  , 
"  njife 

it  follows  that  the  existence  of  a  solution  of  type   (2.52)  implies  the  relation 

(2.61)  sin  J  (a-A)   sin  J(a+A)  D(a,X)  -  0. 
Therefore,   we  find  from    (2.59)  and   (2.60),  by  a  simple  calculation, 

(2.62)  h  sin^  J  A  D(0,X)  -  2  -  2  cos  na  -  2  -  y^(n,X)  -  y2(n,X)   . 

Alternatively,   we  could  have  computed  K  in  terras  of  D(l, X).     The  same  argument  as 
above  would  have  given  us  the  relation 

(2.63)  h  cos^  f  A  D(1,X)  -  2  +  y^in,\)  +  y2(n,X)  . 

Summarizing,  we  have 

Theorem  2.9.     The  discriminant  /\(X)  of  Hill's  equation  (2.1)  can  be 
expressed  in  two  ways  as  an  infinite  determinant  involving  the  Fourier 
coefficients  g     of  Q(x),    (which  are  normalized  so  that  g     -  0  arid 
g_^  ■  g^)}   namely,    with 


D^(X) 


g  r  00 

"n-m  ' ' 


X-Un' 


T  +  6         I 
2         n,m  I 

IP 


and 
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Dj^(X)  • 


g 


n-E 


X-(2n+l) 


?     n.m 


oo 


-00 


We  have: 


In  the  case  where  Q(x)  -  Q(-x),  the  determinants  D  and  D^  can  be  factored  into 
the  product  of  two  infinite  determinants,  each  of  which  can  be  expressed  in  terms 
of  the  factors  y^(|),  ...,  72(5-)  of  ^ -2  and  ^  +  2  (see  Section  1.3,  Theorem  1.3 
and  equation  (2.17)).  We  have 

Theorem  2.10.  Let  e  «  2  for  n  »  1,2,3,  ♦»«  and  e^  -  0.  L»t  sgn  n  -  ♦  1 
for  n  »  1,2,3,  ..»,  sgn  (-n)  ■  -  sgn  n  and  sgn  0-0.  Then  the  four 


infinite  determinants 


C,(X) 


s„(x) 


n.m 


/^T  (X  -  Un  ) 


00 


.5 


+  ^n-ra       ^n-w. 


X  -  Un' 


T 


00 


c^(x) 


n.m 


(gn-m  ^  gn^^l^   I^^^^"  ^  sgn(m^lT) 
/ri-     [\-(2n*l)5 


00 


Si(X) 


^n-m  "  ^n-Hn+l 


00 


X-(2n+l)' 
(where  the  g^  are  real  and  g     -  g^^;  g^  "  O)  satisfy  the  relations 

C^(X)  S^(X)  -  D^(X) 
Cj^(X)  Sj^(X)  -  D^(X) 

Asm  (^A)c^(x)  "y[i^,-k) 
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sin  i  A 

—^  s^(x)  -  y,q.x) 

cos  (n  y/T)C^(X)  -  y^{^,\) 

cos  (n  fk)S^{\)  '   ygCpX) 

For  a  proof  of  Theorem  2.10  see  Magnus,  1955* 

Theorems  2.9  and  2.10  are  useful  for  the  computation  of  the  first  characteristic 
values  of  Hill's  equation.  Theorem  2.9  could  also  be  used  for  the  computation 
of  the  first  terms  ^   (X)  (see  2.U$)  in  the  expansion  of  /\M»     However,  the 
higher  term  will  then  appear  in  a  different  form.  For  example,  we  find  from 
Theorem  2.9  that 

n,a«-oo  (X-n  )  (X-m  ) 

and  it  requires  a  considerable  number  of  calculations  to  derive  from  this  the  result 
in  Corollary  2.6. 

2»U»  Asymptotic  behavior  of t he  characteristic  values 

In  this  section  we  shall  be  concez*ned  with  estimates  (upper  and  lower  bounds) 

for  the  characteristic  values  X  (n  =  0,1,  ...)  and  Xi(ni  -  1,2,3,  ...)  or,  alternatively, 

n  JB 

with  the  location  of  the  intervals  of  stability  as  defined  in  Theorem  2.1  of 
Section  2.1. 

For  large  n  or  m,  Theorem  2.5  and  Corollary  2.6  yield  seme  information  about 
the  X^  and  X^  which,  according  to  Floquet's  Theorem,  are  the  roots  of  the  equations 
/\(X)  -  2  "  0  and  ^(X)  +  2  ■  0.  However,  much  better  results  are  available, 
although  these  are  rather  difficult  to  obtain.  Here  we  shall  merely  state  without 
proofs  a  few  results,  (por  a  detailed  account  of  the  results  available  and  for 
references  to  the  extensive  literature  we  refer  the  reader  to  Starzinskii,  1955 
and  Krein,  1955.)  We  have,  accarding  to  Borg,  19ii)4: 

Theorem  2.11.  Let  Q(x)  in  (2.1)  be  normalized  so  that 


and  let 


i|   |Q(x)| 


-3U- 
dx  -  A. 


o 

Let  n  be  an  integer  such  that 

^  A  <  n  . 

Then  the  (n+l)-st  interval  of  stability  for  (2.1)  contains  the  interval 
defined  by 

or  equivalently  (for  n  ■  1,2,3,  o»»): 

A  result  much  stronger  than  Theorem  2.11  is  known  in  the  symmetric  case  (see 
Section  1.3),  where  the  characteristic  values  are  the  eigenvalues  of  an  ordinary 
boundary  value  problem  of  the  Sturm-Lion ville  type.  We  have  (Borg,  191*6): 

Theorem  2.12.  Let  Q(x)  -  Q(-x),  Q(x+n)  »  Q(x),  and  assume  that  Q(x)  has 
continuous  first  and  second  derivatives.  Let  >._^,  X  ,  X-,  •••,  X-,  X^,..* 
be  defined  as  in  Theorem  2.1.  Let  n  denote  a  positive  integer.  Then  for 


n  — >oo 


X'  -  -  (2n-l)^  +  -^-  +  o(n"2) 
2^^  (Un)2 

x!^  -  (2n-l)2  +  -2—  +  o(n-2) 
^"  (Un)2 


X_  -  Un2  +  _2_-  +  o(n-2) 
2"        (Unf 


X-  ,  -  (Un)2  +  -°   +  o(n-2) 


where 
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2 
C  -  i  I   |Q(x)|  dx 


Theorem  2.12  shows  that,  in  the  symmetric  case  where  Q(x)  ■  Q(-x),  the 
intervals  of  instability  tend  to  zero  like  o(X  )  as  X  — >  co .  No  equivalent 
theorem  seans  to  ba  known  in  the  case  of  a  general  Q(x).  In  fact,  a  result  due 
to  Borg,  19U6,  makes  it  appear  unlikely  for  such  a  theorem  to  hold.  Borg's  Theorem 
states  that  two  sequences  X  and  X*  (n  ■  1,2,3,  •••,  excluding  X  )  determine 
uniquely  an  even  function  Q(x)  such  that  the  X  ,  X«  are  characteristic  values 
for  (2.1)  with  this  particular  Q(x),  provided  that  the  X  and  X'  satisfy  (2.U)    * 
and  certain  asymptotic  conditions  as  n  ->  oo .  If  the  characteristic  values 
belonging  to  (2.1)  for  any  Q(x)  would  satisfy  the  same  asymptotic  conditions  as 
those  belonging  to  an  even  Q(x),  then  it  would  be  possible  to  assign  an  even 
function  Q  (x)  to  every  Q(x)  such  that  (2.1)  would  have  the  same  characteristic 
values  for  either  Q  or  Q.  It  seems  very  unlikely  that  assertion  of  this  type 
should  be  true. 

Theorems  2*11  and  2»12  are  due  to  Borg,  19UU,  19U6.  For  various  related 
theorems  see  Starzinskii,  1955  and  Putnam,  195U.  Of  these  we  shall  formulate 
here  only  the  following  result  which  seems  to  have  been  the  first  of  its  type  and 
which  is  due  to  Liapounoff,  1907: 

Theorem  2.13  Let  p(x)  ^  0  be  a  non-negative  piecewise  continuous 
periodic  function  with  period  n.  Then  all  solutions  of 


y"  ♦  p(x)  y  »  0 
are  bounded  for  all  value a  of  x  if 


p(x)  dx  <  U. 
o 

This  condition  is  best  possible  in  the  sense  that,   for  any  e  >  0,   there 
exists  a  non-negative  piecewise  continuous  function  p   (x)  satisfying 

Pq  3^  0     ,     PQ(x+n)  =  p^(x) 


and 


"I 

f  0 


p^(x)  dx  <  U  +  e  , 
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such  that  at  least  one  solution  of 

y"  +  Pq(x)  y  =  0 

is  unbounded  as  x  -*>■  +  oo  ♦ 

Generalizations  of  Theorem  2.13  which  characterize  the  n-th  interval  of 
stability  of   (2.1)  and  are  best  possible  results  in  the  same  sense  as  Theorem 
2  •13  are  due  to  Borg,   19Uii. 

2«5«     Properties  of  the  solutions 

In  this  and  in  the  following  section,  we  shall  review  some  non-elementary 
results.     Terms  belonging  to  the  general  theory  of  linear  differential  equations 
will  be  used  without  further  explanation. 

The  following  result  has  been  proved  by  Haupt,   191U: 

Theorem  2«lU.     Let  y(x,X)  be  a  non-trivial,   real  periodic  solution  of   (2.1) 
with  period  n  or  2ii.     If  X  -  "K^^^-^  ££  ^  *  ^2n'   ^^^"  ^  has  exactly  2n+l  zeros 
in  the  half  open  interval  0  <  x  <  2n.     If  X  -  X-     ^  or  X  ■  X^  ,   then 
y  has  exactly  2n  zeros  in  0  <  x  <  n. 

If  X  belongs  to  an  interval  of  stability,   any  solution  of   (2.1)  has  infinitely 
many  zeros.     This  conclusion  can  be  derived  from  Floquet's  Theorem  and  from  the 
fact  that  the  general  solution  of   (2.1)  can  be  written  in  the  form 

X  -1 


(2.65) 


/y'  ^72         COS       f     Llit)  *  y2(t)l 


dt  +  a 


where  A,   a  are  constants  and  y.,   y^  are  the  normalized  solutions  introduced  in 
Section  1«2.     Yelchin,   19U6,    showed  how  one  can  decide  whether  a  solution  of  an 
equation  of  Hill's  type  oscillates  or  not.     For  results  and  references  see 
Starzinskii,   1955. 

Periodic  solutions  of  period  it  belonging   to  different  characteristic  values 
X  ,   X     of    (2.1)  are  orthogonal  in  the  interval   (0, n).     Any  two  periodic  solutions 
of  period  n  or  2n  are  orthogonal  in  the  interval   (0,2n)  if  they  belong  to  different 
characteristic  values. 

The  following  result  is  the  analog  of  the  well  known  theorem  about  the  expan- 
sion of  a  function  in  a  Fourier  series.      (See  Weyl,   1910  or  Coddj.ngton  and 
Levinson,   1955).     We  have 
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Theorem  2«l5«     Let  [i  ,   m  »  0,1,2,    ,.«  be  the  characteristic  values 


X  ,   X^,   ^5>^»    •••  in  their     natural  order  and  let  z   (x)  be  an  orthon- 
otmal     set  of  periodic  solutions  of   (2»l)  such  that  z     satisfies  the 
equation 

Then  every  continuous  periodic  function  of  period  2n  whose  second  derivative 
Is  square  Integrable  in  every  finite  interval  can  be  explained  in  a  uniformly 
and  absolutely  convergent  series 

00 

T~    c     z   (x) 

in*C 

with  constant  coefficients  c   . 


The  theory  of  Weyl,   1910,   can  also  be  applied  to  the  differentia]   equation 
(2.1)  for  the  Interval  -oo  <  x  <oo.     There  always  exists  a  solution  of   (2.1)  which 
Is  not  square  integrable  in  (-oo,   o).     An  analogous  statement  is  valid  for  the 
Interval   (o,oo).     We  thus  have  the  limit  point  case  at  both  end  points  -oo  and 
+00.     Since  every  non-trivial  solution  of   (2»1)  fails  to  be  square  integrable  in 
(-co, oo)   (cf.  Floquet's  theorem),   the  spectrum  is  purely  continuous  and  can  be  seen 
to  coincide  with  the  union  of  the  intervals  of  stability  (see  Hartman  and  Wintner, 
19l49)»     Thus,   according  to  Weyl,   1910,  we  have  the  following 

Theorem  2.16.     Let  f (x)  be  a  continuous  and  two  times  differentlable 
function  of  x  which  is  defined  for  -oo  <  x  <oo  and  for  which 


oo  +  oo 


f        if(x)|^dx<oo     ,    I        |f"  (x)|2  dx<oo 

/  -  oo  /-  00 


Let  S  be  the  set  on  the  real  X-axis  which  consists  of  the  union  of  the 


open  intervals  of  stability  of   (2.1)   and  their  endpoints.     Let  y,(x,X) 
and  y-(x, X)  be  the  normalized  solutions  of   (2.1).     Then  there  exist 
functions  M. (X)  and  Mp(X)  defined  on  S  and  such  that 


f  (x)  -  j     .  y^(x,X)  dMj^(X)   +  y2(x,X)  dM2(X) 


According  to  Borg,  19U6,  the  only  case  in  which  the  set  S  is  connected  occurs 
if  Q(x)  £  0.  This  case  leads  to  the  ordinary  Fourier  theorem. 
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The  method  used  in  proving  Theoran  2.2  can  also  be  used  to  show  that  for  every 
fixed  value  of  x  a  non-trivial  solution  y(x,X)   of   (2.1)  is  of  order  of  growth  ^ 
with  respect  to  X.     From  the  estimates   (2.32)  and  from  the  theorem  due  to  Paley 
and  Wiener,   193U   (Theorem  x,    p.   13)  the  following  fact  can  be  derived:     Let  x  be 
real  aiid  let  y(x,X)  be  a  solution  of   (2.1)  for  which 

y(0,X)  -  a     ,     y'(0,X)  -  b     . 

Then  there  exists  a  function  G(x,9)  of  the  real  variables  x,   Q  which  is  defined 
for   |9J   <  |x|   and  is  such  that     for  all  real  values  of  x  and  for  X  >  0 

X 

y(x,X)  -  a  cos(x  A  )  +   [    G(x,e)  e^®  *^  d©       . 

X 

The  function  G-(x,©)  satisfies  the  partial  differential  equation 

^    -  1°     +  Q(x)  G  .  0 

For  details  see  Magnus,   1955,   and  for  applications,   see  Gelfand  and  Levitan,   1955, 
(p.  296). 

2.6.     The  coe?d.stence  probl em 

In  the  «ase  of  an  ordinary  homogeneous  boundary  value  problem  for  a  linear 
differential  equation 

y"  >  ]\*Q{xJ}  y  -  0  , 

(where  Q(x)  need  not  be  periodic)  the  characteristic  values  (eigenvalues)  of  X 
are  all  simple,  which  means  that  there  cannot  exist  two  linearly  independent 
solutions  satisfying  the  same  homogeneous  boundary  conditions  for  a  finite  interval. 
However,  the  condition  that  a  non-trivial  solution  should  be  periodic  (of  period 
n  or  2ft  if  Q  is  of  period  n)  is  of  a  different  nature.  Thus,  we  may  have  double 
characteristic  values  and  two  linearly  independent  solutions  of  period  n  or  2n 
associated  with  the  same  eigenvalue.  The  conditions  for  this  to  happen  have  been 
discussed  in  Section  2.1.  They  can  be  stated  by  saying  that  the  equation 
/\     (X)  -  U  -  0  must  have  a  double  root  or,  more  vaguely,  by  asserting  that  an 

interval  of  instability  must  disappear.  In  this  case  there  exists  an  interval 

2 

of  stability  which  contains  a  double  root  of^  (X)-l4-0o  It  should  be 

emphasized  that  there  always  exist  infinitely  many  characteristic  values 
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X  ,  X,,  X^,  •••  and  X.,  Xp,  X^,  ••.  and  that  th«  set  of  these  contains  the  set 
of  boundary  points  of  the  intervals  cf  instability,  but  that  (with  the  exception 
of  X  )  any  number  of  characteristic  values  may  lie  within  intervals  of  stability. 
The  question  of  the  conditions  which  Q  must  satisfy  to  insure  the  existence  of 
two  linearly  independent  solutions  for  a  certain  set  of  characteristic  values 
shall  be  called  the  coexistence  problem.  The  following  two  important  results 
bearing  on  this  problem  are  due  to  G.  Borg,  19U6: 

Theorem  2.17«  All  of  the  roots  of  the  equation 
^(X)  +  2  »  0 

are  double  roots  (and  two  linearly  independent  solutions  of  (2.1)  with 
period  2n  exist  whenever  one  such  solution  exists)  if  and  only  if  Q(x) 
has  period  n/2. 

Corollary  2.17.  All  of  the  roots  of  the  equation 

^^  (X)  -  U  -  0 

with  the  exception  of  X  the  smallest  one,  are  double  roots  (and  no 
interval  of  instability  exists  for  X  >  X  )  if  and  only  if  Q(x)  is  a 
constant. 

One  half  of  Theorem  2.1?  is  a  consequance  of  the  Corollary  to  Floquet's  Theorem 
in  Section  1,2»     In  fact,   if  Q(x)  is  of  period  n/2,   then  2n  equals  four  times  the 
smallest  period  of  Q(x)  and  therefore  all  of  the  non-trivial  solutions  of   (2,1)  are 
of  period  2n  if  one  of  them  has  this  property.     It  is  much  more  difficult  to  prove 
the  conversp  of  this  statement,   i.e.,   to  prove  the  other  half  of  Theorem  2.17, 
and  we  must  refer  the  reader  to  Borg's  paper  for  the  proof  of  this  result. 

Corollary  2.17  follows  from  Thecran  2.17  i/  we  note  that  Q(x)  would  be  periodic 
with  periods  n/2,   n/h,   n/8,    •••  and  therefore,   a  constant. 

Examples  for  equations  of  Hill's  type  with  a  finite  number  of  intervals  of 
stability  will  be  given  in  the  section  on  Lame^^ (Ince's)  equation. 
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